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Alternate  Approaches  to  Laser  Guided  Discharges 

J.  Davis,  R.  E.  Terry,  G.  Petrov,  A.  Velikovich 
Plasma  Physics  Division,  Naval  Research  Laboratory 
Washington,  D.  C.  20375 


1  Introduction 

The  propagation  of  intense  ultrashort  laser  pulses  in  the  atmosphere  is  currently  under  inten¬ 
sive  investigation  because  of  its  potential  for  a  number  of  applications,  particularly  for  the 
generation  of  super  continuum  or  white  light  radiation,  LIDAR,  sensor  and  detector  areas 
and  as  a  trigger  for  streamers  in  air  for  initiating  lightning  [26], [23], [33], [10], [19].  jn  addition, 
it  has  the  potential  for  directed  energy  and  possibly  for  1ED  neutralization  applications.  A 
desired  property  of  ultrashort  laser  pulses  injected  into  the  air  is  their  potential  to  propagate 
large  distances,  producing  fairly  uniform  ionized  channels  on  their  wake.  This  feature  of 
producing  long  narrow  plasma  channels  is  an  important  property  for  weapons  applications. 
However,  although  the  laser  is  excellent  for  causing  air  breakdown  and  creating  the  plasma 
channel,  very  little  energy  is  transported  to  prospective  targets  to  induce  significant  dam¬ 
age.  A  way  to  improve  this  deficiency  is  to  launch  an  electrical  discharge  [32], [14], [18], [40], [37] 
through  the  plasma  channel.  Laser  guided  discharges  are  under  intensive  investigation  for 
triggering  lightning  in  order  to  divert  the  lightning  from  striking  power  stations.  These  ap¬ 
plications  have  been  discussed  in  a  number  of  papers.  Our  interest  in  laser  guided  discharges 
is  the  potential  to  control  and  to  stabilize  the  directionality  of  the  discharge  path,  while 
maintaining  the  channel  and  projecting  enough  energy  to  cause  damage  to  a  test  object. 
The  concept  that  we  propose  is  based  on  a  complete  circuit,  equivalent  to  a  transmission 
line,  between  a  laser  source  and  a  distant  object.  In  order  to  explore  further  this  aspect  of 
laser  discharge  plasma  channels  it  is  necessary  to  determine  both  the  nature  of  the  interac¬ 
tion  of  the  laser  with  the  atmosphere  and  the  ability  to  deliver  interesting  amounts  of  energy 
to  the  object. 

The  interaction  of  the  laser  with  the  ambient  air  is  much  more  complicated  than  originally 
supposed^-*-]’^’^^’^^’^®]’^^.  As  the  laser  interacts  with  and  penetrates  the  atmosphere 
it  can  filament  into  a  number  of  streamers  reducing  the  power  delivered  to  a  distant  object. 
We  propose  to  model  and  numerically  simulate  the  atmospheric  interaction  and  response  of 
the  atmosphere  to  the  propagation  of  an  intense  ultrashort  IR  laser  pulse.  The  atmospheric 
response  to  the  pulse  will  be  modeled  self-consistently  coupled  to  the  equations  representing 
the  chemistry  of  reaction  kinetics  of  the  atmospheric  constituents  as  the  quiescent  air  is 
ionized  and  follow  its  subsequent  return  to  ambient  conditions.  It  is  crucial  to  determine  the 
lifetime  and  the  length  of  plasma  channels  as  they  develop  and  evolve  if  a  pulsed  IR  laser 
produced  plasma  channel  is  to  persist  long  enough  to  launch  an  electrical  discharge  through 
the  plasma  channel.  The  modeling  involves  numerically  solving  the  propagation  equations 
including  a  number  of  important  processes  such  as  the  nonlinear  Kerr  effect,  self-focusing, 
Raman  scattering,  multiphoton  /  tunneling  ionization,  air  ionization  (as  an  energy  loss)  and, 
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in  the  case  of  laser  triggered  discharges,  air  chemistry.  The  equations  of  propagation  have 
been  solved  for  a  number  of  initial  conditions  and  are  well  documented  in  the  literature. 
However,  for  much  of  this  research  a  variety  of  approximations  are  introduced  to  avoid  the 
complicated  interaction  between  the  laser  deposition,  hydrodynamic  response^,  and  the 
complex  of  reactions  occurring  in  the  air.  Real  air  contains  aerosols  and  water  vapor  ^  that 
can  severely  affect  propagation  lengths.  They  can  act  as  nuclcation  points  for  ionization  and 
create  strong  plasma  inhomogeneities.  Aerosols,  dust  and  other  warfighter  environments 
therefore  pose  additional  problems  causing  laser  energy  absorption  and  scattering  and  al¬ 
lowing  air  breakdown  to  occur  at  even  lower  laser  intensities. 

The  concept  advanced  here  offers  a  coordinated  path  to  the  practical  exploitation  of  laser 
guided  discharges  as  energy  delivery  systems.  We  discuss  the  simplest  designs  for  systems 
of  this  type  as  well  as  the  modeling  efforts  and  experimental  programs  required  to  reduce 
them  to  practice. 


2  Background 


In  order  to  better  understand  the  first,  high  voltage  phase  of  the  concept  and  its  relation  to 
lightning,  the  mechanism  of  a  lightning  flash  is  briefly  summarized  [30],[12] _ 

Clouds  can  store  a  tremendous  amount  of  electrostatic  energy,  which,  under  the  right  circum¬ 
stances,  can  be  released  in  a  lightning  flash.  The  lightning  flash  is  initiated  by  an  electrical 
breakdown  process  in  the  cloud,  called  preliminary  breakdown.  This  process  leads  to  the 
creation  of  a  column  of  charge  called  stepped  leader  that  travels  from  cloud  to  ground  (Fig. 
1).  On  its  way  to  the  ground  the  leader  may  give  rise  to  several  branches. 

As  the  stepped  leader  approaches  the  ground  the  electric  held  at  ground  level  increases. 
When  the  stepped  leader  gets  close  (a  few  hundred  meters)  the  electric  held  at  the  grounded 
structure  increases  to  such  a  level  that  electrical  discharge  is  initiated  from  them.  This 
discharge,  called  connecting  leader,  travel  toward  the  down-coming  stepped  leader  (Fig.  2). 
The  connecting  leader  may  successfully  bridge  the  gap  between  the  ground  and  the  down¬ 
coming  stepped  leader  (Fig.3).  The  object  that  initiated  the  connecting  leader  is  struck  by 
lightning. 

Once  the  connection  is  made  between  the  stepped  leader  and  ground,  a  wave  of  near  ground 
potential  travels  along  the  channel  toward  the  cloud.  This  electromagnetic  wave  travels  at  a 
velocity  of  nearly  the  speed  of  light.  This  is  called  the  return  stroke  (Fig. 4)  and  it  causes  the 
actual  damage.  The  diameter  of  the  lightning  return  stroke  is  typically  several  centimeters. 

The  available  data  indicates  that  the  peak  current  is  10-100  kA  and  the  energy  delivered  to 
the  stricken  object  vary  from  several  MJ  to  about  1  GJ. 
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100MV 


Figure  1:  Downward  moving  step  leader  travels  from  cloud  to  ground 
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Figure  2:  A  connecting  leader  moving  upward  is  launched  from  a  tall  object. 
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Figure  3:  The  connecting  leader  meets  the  stepped  leader  from  the  cloud. 


Figure  4:  The  return  stroke. 
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Since  the  1970s,  many  researchers  have  studied  the  triggering  and  guiding  of  electrical  dis¬ 
charges  in  laboratory  with  different  types  of  laser  systems.  The  first  investigations  of  laser 
control  of  large-scale  discharges  have  been  done  primarily  with  long  pulse  IR  lasers.  Nev¬ 
ertheless,  the  use  of  such  lasers  suffers  from  important  limitations,  leading  researchers  to 
consider  other  candidates,  such  as  long-pulse  ultraviolet  lasers  or,  more  recently,  ultrashort 
pulse  lasers.  The  advantage  of  ultrashort  laser  pulses  is  their  very  high  intensity,  which  al¬ 
lows  the  creation  of  long  continuous  plasma  channels  in  air  via  multiphoton  ionization  with 
a  modest  amount  of  energy.  The  plasma  channel  is  used  to  propagate  and  guide  a  plasma 
wave  (leader).  The  latter  is  a  highly  conductive  media,  which  travels  from  one  electrode 
toward  the  other.  Once  the  leader  bridges  the  gap  between  the  electrodes,  a  wave  of  poten¬ 
tial  travels  along  the  channel  with  a  speed  close  to  the  speed  of  light  delivering  MJ  levels  of 
energy. 


3  Theoretical  Description  and  Operation 


The  damaging  behavior  of  natural  lightning  can  be  enhanced  by  adding  a  final  pulse  power 
element:  closing  a  circuit  and  delivering  a  tailored  pulse,  similar  to  the  action  of  a  ”TASER” 
device.  The  reason  for  having  a  four-step  process  (laser  created  channel-leader-return  stroke- 
final  pulse)  instead  of  just  one  is  by  far  dominated  by  one  factor,  viz.  the  conductivity  of 
the  media.  Only  a  medium  with  very  high  conductivity  is  good  for  delivering  a  lethal 
amount  of  energy  to  the  target.  The  plasma  filaments  created  by  the  laser  (step  one)  have 
low  conductivity.  During  step  two  (leader  propagation)  the  conductivity  of  the  channel 
increases  by  several  orders  of  magnitude.  During  step  three  (the  return  stroke),  the  channel 
conductivity  increases  even  further,  allowing  swift  and  efficient  delivery  of  energy  to  the 
target.  During  step  four  (pulse  delivery)  the  ohmic  dissipation  along  the  line  drives  the 
conductivity  up  even  further.  The  individual  steps  will  be  outlined  next. 


3.1  Channel  Formation 

The  first  part  of  the  problem  is  to  create  a  plasma  channel  ^102-103  m  long.  Ultrashort  pulse 
lasers  (USPL)  are  the  best  candidates  for  the  creation  of  long  continuous  plasma  channels 
in  air  Fl>[6], [5], [25], [13], [36], [4], [8], [28], [22], [16], [3].  Nonlinear  propagation  of  ultrashort  pulses  in 
air  efficiently  creates  free  electrons  through  multi-photon  ionization  of  oxygen  and  nitrogen. 
A  plasma  channel  («0.2  mm  in  diameter)  is  formed  with  a  typical  electron  density  of  ne 
RilO15  -1018  cm-3  and  electron  temperature  of  1-10  eV.  This  plasma  channel  can  be  very 
long,  hundreds  of  meters.  For  example,  the  generation  of  long,  stable  plasma  filaments  in  air 
over  distances  in  excess  of  200  m  has  been  reported^.  The  possibility  of  creating  plasma 
channels  hundreds  of  meters  long  is  critical  for  the  present  application  since  these  filaments 
can  provide  a  preferred  channel  for  long  streamer  discharges  and  guide  leaders. 
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The  energy  requirements  can  be  estimated  as  follows.  Most  of  the  laser  energy  is  spent  for 
optical  field  ionization  (OFI), 


&OFI  neVEion  ,  (1) 

where  V  =  ncPL/A  is  the  plasma  volume  and  Eion  is  the  ionization  energy.  For  typical  values 
nepsl016cm~3,  d  =  200  /zrn,  L  =  100  -  1000  m  ,  and  Eion  «  1.2eV  =  1.6  •  10-18  J  the  laser 
energy  spent  for  optical  field  ionization  is 


Sofi  ~  0.06  -  0.6[J]  .  (2) 

These  energy  requirements  are  not  too  restrictive,  but  are  certainly  still  near  the  upper  limit 
of  the  presently  available  USPL  systems.  What  is  of  primary  interest  is  the  time  evolution  of 
the  channel.  After  only  «  0.1  /zs  following  the  plasma  channel  formation,  the  electron  density 
(and  conductivity)  decreases  to  ne«1013cm-3  and  keeps  decreasing  in  time.  The  electron 
population  doesn’t  really  disappear,  but  attaches  to  oxygen,  so  that  an  ion  population 
arises,  which  figures  heavily  in  the  physics  of  leader  propagation  and  the  subsequent  return 
stroke  -  a  generic  population  history  of  these  species  is  shown  in  Figure  5.  According  to 
some  analyses,  the  plasma  channel  retains  properties  which  are  still  adequate  for  triggering 
and  guiding  of  the  discharge  for  tens  of  /zs,  but  ultimately,  the  channel  will  lose  its  guiding 
ability. 


Figure  5:  Conductors  in  the  Laser  Induced  Breakdown  channel. 


3.1.1  Propagation 


The  model  for  long-distance  propagation  is  a  2D  version  of  the  nonlinear  Schrodinger  equa¬ 
tion  for  the  envelope  of  the  electric  field.  Our  starting  point  is  the  wave  equation  for  the 
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electric  field  E 

(vi  +  ^-^)E  =  S,  (3) 

where  is  the  transverse  Laplacian,  c  is  the  speed  of  light,  z  is  the  direction  of  propagation 
and  S  is  a  source  term.  The  laser  electric  held,  as  well  as  the  source  term  are  written  as 
complex  amplitudes  multiplied  by  a  rapidly  varying  phase  -0  =  k0z  —  u0t\ 

E(r,  z,  t )  =  A(r,  z,  f)e!^’^ex/ 2  +  c.c.  (4) 

S (r,z,t)  =  S(r,  z,  f)e*^z,^ex/2  +  c.c.  (5) 

Here  k0  =  ujq/c  =  2tt/X  is  the  wave  number,  cc0  is  the  laser  frequency,  ex  is  the  unit  vector 
in  the  direction  of  polarization  and  c.c.  denotes  complex  conjugate. 

The  source  term  in  the  right  hand  side  of  Eq.  (3)  consists  of  several  components,  which 
account  for  the  action  of  the  media  through  which  the  laser  propagates.  For  the  laser 
intensity  and  laser  pulse  duration  of  interest,  the  target  media  can  be  described  as  partially 
ionized  plasma.  More  specifically,  it  consists  of  nitrogen  and  oxygen  molecules  (the  air 
constituents),  their  dissociation  products,  and  the  plasma  that  is  generated  by  the  leading 
edge  of  the  laser.  For  N2-O2  mixture  at  atmospheric  pressure  and  a  typical  ionization  degree 
of  ICE'  —  10-3  the  source  terms  read: 


(6) 

(7) 

(8) 

(9) 


(v)  ionization 

TTi°n  Q  Tjion  Q 

=  -ik"°^A  .  (10) 

c\A\z  or  I  or 

There  are  several  parameters  participating  in  the  source  terms:  n0  is  the  refractive  index 
of  air,  assumed  to  be  unity,  n2  is  the  nonlinear  refractive  index  of  air,  /  =  ^  is  the  peak 
laser  intensity,  Xl  is  the  laser  susceptibility  at  frequency  uq,  Q  is  the  response  function  for 
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stimulated  Raman  scattering,  tir  is  an  effective  nonlinear  refractive  index  due  to  Raman 
scattering,  v  is  electron- neutral  collision  frequency,  kp  =  ujp/c  is  the  plasma  wave  number 
with  plasma  frequency  up  =  (47T e2ne/m)1/2,  e  and  m  are  the  electron  charge  and  mass 
respectively,  ne  is  the  electron  density,  e\A\/(mcuJo)  is  the  normalized  electric  held  strength, 
Ulon  is  the  ionization  potential  of  O2,  and  dne/dr  is  the  OFI  rate. 


An  equation  suitable  for  numerical  solution  is  obtained  after  a  sequence  of  intermediate 
steps.  First,  the  electric  held  (4)  and  source  term  (5)  are  inserted  into  the  wave  equation 
(3)  and  the  rapidly  varying  exponential  factor  is  canceled  from  both  sides  of  the  equation. 
The  resulting  equation  is  further  transformed  into  the  so  called  laser  frame  by  substituting 
the  variables  (z,t)  with  variables  (z,r),  where  r  —  t  —  z/c  is  the  retarded  time.  Finally,  the 
equation  is  written  in  the  form 


d  A 

Vj A  +  2iko— — 
oz 


d2A  d2A  2,  a 
+  wx  —  c  kofo 


dzd(cr )  dz 2 


d2A 

3{ct)'a 


+  MA  =  0  , 


(11) 


where 

M  =  fc»("2 + ”s)/  _  k2”  0  _  '0) + (2^)  +  ika nrw 

combines  the  non- derivative  terms  of  the  equation,  and  the  transverse  Laplacian  in  cylindrical 
geometry  reads  VyA  =  1J1  (r^0).  The  term  c2fco/32 ofc^)2  ^ie  so_called  Group  Velocity 
Dispersion  (GVD),  which  accounts  for  the  frequency  dependence  of  the  refractive  index  of 
air.  In  the  transformed  coordinates  time  r  refers  to  the  front  of  the  laser  pulse.  Therefore  A 
is  always  bounded  between  0  and  R  in  radial  direction  and  0  and  T  in  the  time  domain,  where 
R  is  an  appropriate  distance  from  the  axis,  much  larger  compared  to  the  laser  pulsewidth, 
and  T  is  the  laser  pulse  duration.  Equation  (11)  is  a  partial  differential  equation,  which  is 
solved  using  finite  differences  by  advancing  the  laser  electric  held  amplitude  A(r,  r)  forward, 
in  +z  direction. 


Initial  and  boundary  conditions  are  required  to  solve  Eq.  (11).  The  initial  condition  at 
z  —  0  comes  from  the  imposed  laser  intensity  profile  /0(r,  r)  upon  entering  the  computational 
domain.  The  input  laser  held  is  then 

A(r,  z  =  0,  t)  =  (8vr/0(r,  r)/c)1/2  ,  (12) 

while  the  boundary  conditions  read: 

d A(r  =  0,  z,  r)/dr  —  0  , 

A{r  =  R,z,t)  =  0  , 

A(r,  z,t  =  0)  =  0  , 

A(r,  z,  r  —  T)  —  0  . 


The  laser  propagation  in  air  creates  plasma  through  OFI.  The  plasma  hlaments  have  local 
electron  density  determined  by  the  electron  continuity  equation 

dne 
dr 


WoFi‘no2  ■ 


(13) 


with  boundary  condition  ne(r,z,T  =  0)  =  0.  In  (13)  the  terms  Wofi  and  no2  are  the  OFI 
rate,  which  is  a  strong  function  of  the  laser  intensity  and  the  density  of  oxygen  molecules, 
respectively.  The  plasma  generation  is  included  into  the  propagation  equation  via  three 
terms,  (8,9,10).  The  real  part  of  (8)  and  the  term  (9)  are  ’’conservative”,  i.e.  they  only 
focus  or  difract  the  laser  beam  without  causing  any  loss  of  energy.  The  imaginary  part 
of  (8),  as  well  as  the  last  term,  (10),  are  depleting  the  laser  beam  energy  by  heating  the 
electrons  through  Inverse  Bremsstrahlung  (IB)  and  ionization,  respectively.  The  energy 
balance  equation,  derived  from  Eq.  (11),  reads: 


dE 

dz 


2 

OT 


o^(r> 

uj20c 


I(r)\drdr. 


(14) 


Equation  (14)  states  that  as  the  laser  advances  it  loses  energy  in  ionization  (the  first  integrand 
term  on  the  right  hand  side  of  (14))  and  IB  (the  second  integrand  term).  Equations  (11)-(13) 
form  a  complete  system  of  equations  for  the  laser  electric  held,  laser  intensity  and  electron 
density. 


3.2  Leader  Propagation 

Leader  propagation  guided  by  a  high  intensity  USPL  has  been  conclusively  demonstrated 
by  several  research  groups.  For  example,  one  of  them  has  performed  an  extensive  study 
of  laser-induced  discharges  in  the  classical  configuration  of  rod-plane  gaps  of  several  me¬ 
ters,  in  appropriate  conditions,  and  obtained  free,  stable  propagation  of  positive  leaders. 
These  experiments  showed  that  a  plasma  channel  created  with  an  USPL  can  trigger  leader 
propagation  at  a  reduced  voltage  and  guide  it  on  a  distance  of  a  few  meters  [20], [34], [29], [2] _ 

The  leader  is  a  thin  («mm)  highly  conductive  filament  that  starts  from  the  active  electrode 
and  propagates  toward  the  grounded  electrode  (or  earth).  The  strong  electric  held  in  front 
of  the  leader  ionizes  the  gas  and  creates  streamers,  which  the  leader  engulfs  and  advances.  In 
general,  leaders  can  be  incepted  and  propagate  in  air  without  a  plasma  channel  (lightning  is 
a  typical  example).  In  this  situation  the  leader  path,  however,  is  to  some  degree  random  and 
the  leader  can  go  virtually  anywhere.  Therefore,  in  order  to  successfully  bridge  (and  short- 
circuit)  the  gap  between  high-voltage  source  and  an  object,  the  leader  must  be  ’’guided”. 
Hence  the  plasma  channel  formation  is  an  integral  part  of  the  process.  There  are  additional 
benefits  of  the  plasma  channel.  The  presence  of  the  channel  alleviates  the  conditions  for 
leader  inception  (lower  breakdown  voltage)  and  increases  the  leader  velocity  by  a  factor  of 
ten  compared  to  natural  leader.  The  first  benefit  directly  affects  the  requirements  and  design 
of  the  power  supply.  The  second  one  may  be  a  lot  more  important  for  long-distance  guiding 
of  leaders  as  the  leader  channel  may  be  capable  of  guiding  for  only  a  limited  period  of  time 
(10  -  100  /j, s)  and  higher  leader  velocity  means  longer  guiding  distance. 

Since  the  whole  process  of  delivering  a  lethal  amount  of  energy  to  the  target  is  driven  by 
the  channel  conductivity,  it  is  instructive  to  compare  the  plasma  conductivity  of  the  laser 
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created  channel  (step  one)  and  the  leader  (step  two).  Such  a  comparison  would  explain  the 
necessity  of  this  stage.  During  step  two  (leader  propagation)  the  channel  resistance  decreases 
by  several  orders  of  magnitude,  as  it  can  be  seen  from  the  following  table. 


Table  1:  Channel  Parameters  Compared 


Parameter 

plasma  channel 

leader 

return 

Length  [m] 

100-1000 

100-1000 

100  -  1000 

Diameter  [cm] 

0.02 

0.2 

2.0 

Gas  temperature  [K] 

300 

3000 

20000 

Neutral  density  [cm-3] 

2.5xl019 

2.5xl018 

5.0xl017 

Electron  density  [cm-3] 

o 

CO 

o 

CO 

5.0xl017 

Conductivity  [fl-1  cm-1] 

O 

1 

CO 

O 

1 

to 

o 

1 

Resistance  [0  cm-1] 

3xl06 

3x10s 

3.0xl0-3 

Current  [A] 

3xl0-3 

3 

3000 

As  can  be  seen  from  the  table,  during  the  second  step  the  channel  resistance  per  unit  length 
decreases  by  a  factor  of  thousand,  both  due  to  increased  conductivity  and  increased  channel 
cross  section  S  =  Tid2 /4.  The  current  flowing  through  the  channel  increases  accordingly, 
which  causes  further  gas  heating  and  channel  expansion.  At  the  end,  the  conductivity 
suitable  for  the  ’’return  stroke”  is  realized. 


3.3  Pulse  Power 

The  energy  delivery  technology  of  interest  comprises  two  subsystems:  a  high  voltage  leader 
propagation  interface  and  a  low  voltage  transmission  line  for  final  pulse  delivery. 


3.3.1  High  Voltage 


The  schematic  representation  of  the  transmission  line  interface  to  the  systems  described 
above  is  shown  in  Figure  6.  The  plane  parallel  capacitor  faces  are  shown  perforated  by  a 
laser  beam  convolute.  A  high  breakdown  strength  dielectric  occupies  the  space  immediately 
adjacent  to  the  convolute  and  is  shown  there  with  a  gas  gap  at  NTP  between  it  and  the  laser 
beam  convolute. 

In  the  early  high  voltage  phase  of  operation  the  initial  laser  channels  are  quite  resistive  (1.2  - 
3.0  kfl/cm)  and  support  only  minimal  currents.  Connected  to  the  hot  side  of  their  respective 
capacitors  through  a  convolute,  an  assessment  must  be  made  of  the  leakage  currents  to  the 
ground  plane  of  the  capacitor  through  the  gas  gap  around  the  laser  beam  convolute.  If  these 
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PLANE  PARALLEL  DISK  LINE  WITH  A  LASER  BEAM  CONVOLUTE 


axial  laser 
channel 


cathode 


Figure  6:  Transmission  Line  Interface 


currents  are  too  great,  leader  propagation  will  not  be  supported  and  a  hard  electrical  contact 
to  the  load  will  fail. 

Making  the  approximation  that  the  laser  beams  don’t  heavily  perturb  the  capacitor  fields 
where  the  convolutes  are,  the  expected  leakage  currents  can  be  estimated  from  those  fields 
and  an  estimate  of  gap  conductivity  alone.  As  shown  in  Figure  7,  which  uses  the  analytic 
unperturbed  electrostatic  potential  solution  (in  planar  oblate  spheroidal  coordinates), 


(C cot  HO-i) 


mo  =  Eo  <[C 


(15) 


7 T 


with  a  characteristic  length  (a,  the  hole  radius)  and  Eq  a  characteristic  field  strength  for  the 
capacitor,  see  also  [31].  Examining  the  field  components,  the  field  E ^  is  the  radial  component 
near  the  axis.  Once  the  laser  beam  convolute  is  made  small  compared  with  the  hole  size, 
this  mostly  radial  field  component  will  be  only  about  0.05  E0.  The  orthogonal  component, 
mostly  axial,  aligns  with  the  unperturbed  field  of  the  plane  parallel  disk  feed  far  from  the 
hole. 

Using  a  typical  electron  density  of  corona  discharge  in  air  at  perhaps  about  1011  cm'3,  the 
expected  resistance  for  the  parasitic  radial  current  can  be  approximated  by 


(16) 
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Figure  7:  £  Field  Strength  Contours 


with  fgeo  a  factor  of  order  one,  and  lgap  the  physical  size  of  the  capacitor  separation.  When 
such  an  estimate  is  used  in  a  parallel  resistive  current  division  of  channel  current  and  parasitic 
current,  the  expected  leakage  current  ratio  can  be  written, 


*  leak 
I  chan 


,0.01 


nP 


Jgeo 


){LgaplCap){^ll  ) 


(17) 


The  held  pattern  just  described  around  the  convolute  can  only  persist  for  the  duration  of 
the  high  voltage  phase.  As  a  return  stroke  along  the  laser  channels  brings  the  conductivity 
up  by  several  orders  of  magnitude,  this  held  pattern  will  shift  to  that  of  a  conventional 
post  hole  convolute  with  a  larger  radial  held  strength.  Balancing  that  radial  electrical  stress 
will  depend  upon:  (i)  the  time  developing  drop  in  the  voltage  across  the  disk  feed,  (ii)  the 
increase  of  magnetic  insulation  around  the  laser  channel,  and  (iii)  the  dielectric  strength  of 
the  feed.  The  end  state  of  that  process  is  a  disk  feed  at  ground  potential  perforated  by  a 
highly  conductive  laser  channel. 


3.3.2  Laser  Channel  Transmission  Line 


Once  formed  the  laser  channels  and  their  connected  circuit  through  a  load  can  be  considered 
as  a  twin  lead  or  Lecher  line  -  described  by  macroscopic  specihc  capacitance,  inductance, 
and  series  or  shunt  resistances  as  shown  in  Figure  8. 
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Owing  to  the  resistance  per  unit  length  p  in  the  laser  channels  forming  this  line,  with  a 
separation  h  and  a  conductor  diameter  d,  it  is  described  by  a  line  impedance 


Z,ine  =  1201n|  t  +  \/(t)2_l  |  [O], 


(18) 


for  a  finite  sized  cylindrical  conductor  pair  and  the  lossy  Telegrapher’s  equation  for  voltages. 

With  CJ  (x,  £)  V forward  T  V backward,  ? 


d2g  (x,  t)  +  ^  dg  ( x ,  t)  _  d2#  (x,  t) 


<9x2 


<9x 


2<9f2 


(19) 


disk 


=  0 


V(t)  I  PARALLEL  LASER  BEAM  TRANSMISSION  LINE 


Figure  8:  Final  Transmission  Line 


The  series  resistance  along  the  plasma  channel  conductors  is  now  represented  by  the  spatial 
damping  decrement,  <5,  which  is  the  channel  resistance  per  unit  length  p  divided  by  the 
line  impedance,  Z\ine.  Explicit  losses  due  to  shunt  resistances  along  the  line,  were  the  two 
channels  to  establish  a  conducting  path  between  them,  would  be  described  by  an  additional 
temporal  damping  term. 

The  forward  and  backward  signals  (c.f.  Figure  9)  on  a  uniform  lossy  line  of  this  nature 
can  be  represented  by  the  following  fundamental  solutions,  parameterized  by  the  free  space 
wavelength  A, 

Vforward  =  COS  ^2  7T  ~  T  ^1  +  1/16  ^  ^  ^  e“1/2r75  (20) 

Vbackward  =  ~  COS  ^2  7T  (^  +  T  ^1  +  1/16  e“1/2?'(“_€)  (21) 
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Kill  Pulse  Waveforms 


Figure  9:  Transmission  Line  Waveforms 


where  £  =  x/X,  r  =  ut,  and  DLoad  represents  the  distance  to  the  load. 

As  shown  in  these  equations,  three  fundamental  dimensionless  parameters  k  =  V  = 

-^7-,  and  5  =  DL°ad  span  the  solution  space.  The  effective  reflection  coefficient  at  the 

terminating  resistive  load  is  just  and  is  invariant  for  impedance  ratios  which  are  the 
inverse  of  one  another.  Everything  of  relevance  for  power  transfer  to  a  load  can  be  specified 
as  a  function  of  [5, 77,  ft],  e.g.  the  coupling  fraction  shown  in  Figure  10  as  function  of  n  with 

[5=10,  77=0.05].  The  line  also  admits  some  dispersion  in  that, 


u(k0)—uj0J  1  +  1/16^ — ^  (22) 

y  ^ line  7T 

as  shown  in  Figure  11.  Under  most  conditions,  this  dispersion  doesn’t  lower  the  group 
velocity  significantly. 


14 
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Figure  10:  Transmission  Line  Coupling 


Dispersion  and  Group  Velocity 


Dispersion  Factor 
Group  Velocity 


Figure  11:  Dispersion  Characteristics 
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4  Proof  of  Principle  Experiments 


An  operational  system  is  thus  decomposed  into  the  four  distinct  phases  discussed  above.  It 
is  certainly  possible  to  test  the  behaviors  of  each  component  in  an  experiment  of  modest 
distance  (~100  m)  and,  once  successful,  the  development  to  a  system  of  interest  to  the 
services  is  then  only  a  factor  ten. 


4.1  Laser  Propagation  and  Channel  Formation 

There  exist  currently  available  laser  systems  which  can  be  used  for  the  initial  experiment. 
One  commonly  known  variety  is  the  Ti-sapphire  system,  which  might  offer  the  parameters: 

•  wavelength  -  795  nm, 

•  energy  -  10  to  50  mJ, 

•  pulse  duration  -  100  to  1000  fs, 

•  power  -  10  to  100  GW, 

•  and  filament  lengths  of  200  m. 


-40  300 


Figure  12:  Initial  Laser  Profile 


The  particular  test  case:  I0  =  1.5  •  1013  W/cm2,  A0=800  nm,  r0=20  /irn,  T=300  fs,  and 
Eo=3  mJ  has  the  initial  profile  shown  in  Fig.  12.  Such  a  pulse  is  advanced  over  space  and 
time  in  2D  via  the  propagation  model  (§3.1.1). 
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Figure  13:  Laser  Profile  at  lm 
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Figure  14:  Laser  Attenuation  with  Distance 
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The  expected  state  of  the  laser  filament  after  a  meter  of  propagation  is  shown  in  Fig.  13, 
while  the  attenuation  of  the  energy  contained  in  the  pulse  is  shown  in  Fig.  14.  The  observed 
damping  decrement  in  this  calculation  is  about  0.32  mJ/m  in  agreement  with  other  work. 
Here  the  model  works  with  only  a  single  filament,  a  very  idealized  propagation  mode.  Yet, 
if  this  damping  were  to  hold  over  10m  lengths,  then  a  50  mJ  input  pulse  would  still  deliver 
18  mJ  in  its  envelope  at  100  m.  Such  an  intensity  would  be  sufficient  to  support  ionization 
throughout  the  channel. 


4.2  Leader  Propagation 

To  establish  a  good  leader  in  a  100  m  channel  to  a  remote  load,  the  high  voltage  pulse 
power  subsystem  will  need  to  erect  a  bipolar  30-50  MV  signal  on  each  laser  channel.  With 
the  typical  resistances  of  the  post  leader  medium  at  about  3.0  kfl/crn  the  early  discharge 
currents  will  be  about  1  -  2  A  and  transfer  only  fraction  of  a  Coulomb  across  the  gap  to  set 
the  channel. 


4.3  Return  Stroke  Connection 

As  the  return  stroke  discharges  the  voltage  brought  across  by  the  leader,  that  connected 
channel  breaks  down  to  a  quite  lower  impedance  of  perhaps  0.3  O/rn  or  less.  Now  for  the 
low  voltage  100  m  pulseline  one  needs  only  perhaps  300  -  500  kV  DC  to  draw  about  10  kA 
over  the  gap  and  help  to  ’’burn  in”  the  channel  even  more.  For  these  low  frequency  signals, 
the  linear  damping  decrement  5  should  be  modest  but  not  severe.  Frequencies  of  about 
300  kHz  should  be  about  right  for  this  component  of  the  signals  sent  up  the  laser  channel 
transmission  line. 


4.4  Pulse  Launch,  Control,  and  Delivery 

For  the  higher  power,  more  damaging,  signals  the  frequency  components  should  be  in  the 
domain  of  30  MHz  or  higher.  For  this  domain  the  damping  decrement  will  be  less  than  or 
about  0.05  and  good  coupling  can  obtain.  The  final  design  of  the  low  voltage  pulseline  driver 
will  therefore  require  some  detailed  attention  to  pulse  shaping. 

The  damage  to  various  loads  can  be  expected  to  depend  upon  several  mechanisms  l^bH Y],[24] _ 
Simple  Joule  heating  of  a  conducting  skin  may  ignite  secondary  energy  release  within  the 
load.  Energy  deposited  in  a  plasma  blanket  evolved  from  a  load  surface  may  develop  a 
significant  impulse  delivery  to  the  load.  Hence  the  experiment  must  examine  energy  delivery 
at  a  short  distance  first,  even  as  small  as  a  meter,  and  start  to  build  a  database  of  material 
responses. 
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5  New  Features  and  Advantages 


The  principal  innovations  we  advocate  are  (i)  the  use  of  a  completed  circuit  of  two  laser 
channels  rather  than  a  single  channel,  (ii)  the  construction  of  both  high  voltage  and  low 
voltage  pulse  power  subsystems  for  channel  breakdown  and  final  pulse  delivery,  (iii)  a  laser 
beam  convolute  to  make  connections  to  both  the  pulse  power  subsystems  within  a  single 
element,  and  (iv)  the  maintenance  of  laser  energies  high  enough  for  optical  field  ionization 
but  lower  than  the  fields  required  for  avalanche  ionization. 

The  dual  channel  arrangement  serves  to  form  a  transmission  line  of  low  and  controlled 
impedance,  in  contrast  to  any  approach  to  sending  power  down  a  single  laser  channel.  With 
only  one  conductor  the  system  inductance  is  markedly  larger  and  the  frequency  window 
for  exciting  low  voltage  pulses  is  thus  curtailed  to  lower  frequencies.  The  laser  twin  lead 
geometry  admits  essentially  all  RF  frequencies  as  TEM  pulses  so  that  signals  tailored  to  the 
target  can  be  constructed  if  desired. 

By  convoluting  separate  high  and  low  voltage  pulse  power  systems  into  the  laser  beams  we 
can  solve  two  easier  electrical  engineering  tasks  instead  on  one  difficult  task.  Each  pulse 
system  can  be  configured  to  operate  on  independent  time  bases  and  voltage  domains  with 
conventional  technology. 

Finally,  by  relying  on  optical  field  ionization  by  fast  pulse  lasers  and  feeding  the  laser  energy 
into  the  channel  in  a  measured  way,  we  create  a  cold  gas  channel  that  will  not  tend  to  bloom 
out  the  laser  beams.  Moreover,  by  using  the  long  lived  species  O7  as  the  charge  reservoir, 
we  can  gain  access  to  a  physics  domain  capable  of  lowering  the  breakdown  voltages  needed 
to  strike  a  low  impedance  lightening- like  arc  to  the  load. 
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